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We suggest a way of rationalizing an intra-seasonal oscillations (IOs) of the Earth atmospheric 
flow as four meteorological relevant triads of interacting planetary waves, isolated from the system 
of all the rest planetary waves. Our model is independent of the topography (mountains, etc.) and 
gives a natural explanation of IOs both in the North and South Hemispheres. Spherical planetary 
waves are an example of a wave mesoscopic system obeying discrete resonances that also appears 
in other areas of physics. 



Introduction. Concept of mesoscopic systems most of- 
ten appears in condensed matter physics, e.g, in studying 
properties of superconductors on a scale comparable with 
that of the Cooper pairs of miniaturized transistors 
on a computer chip, of disordered (glassy, granular) sys- 
tems, when self-averaging is inefficient and fluctuations 
or the system prehistory become important. Similar sit- 
uation occurs also in various natural phenomena - from 
wave turbulent systems in the Ocean [2j and Atmosphere, 
when wave lengthes are compatible with the Earth ra- 
dius 0], to medicine and even in sociology and eco- 
nomics, when finite size of a system (population, sociolog- 
ical group, market) becomes important Mesoscopic 
regimes are at the frontier between detailed, dynamical 
and statistical, self-averaging description of systems. Im- 
portant observation for finite-size, weakly-decaying wave 
systems was made in |: discrete spacial-time resonances 
form small isolated clusters of interacting modes without 
energy exchange between the clusters. Clearly, there is 
exists relatively short, "threshold" wave, involving into 
the cluster with size large enough to "penetrate" into re- 
gion of very short waves, where statistical description, 
that ignores the resonance discreetness, is valid. 

Mathematical problem of finding these clusters in con- 
crete cases is equivalent to solving some systems of high 
(> 12) order Diophantine equations on a space of 6-8 
variables in big integers Vft . Recently developed algo- 
rithms for their analysis [8[ allows to find, in particular, 
all resonance clusters of atmospheric planetary waves, de- 
scribed by the spherical functions Y™ with eigenvalues 
|m| < I < 50. It turned out that in this domain, consist- 
ing of 2500 spherical eigen-modes Y™, there exist only 
20 different clusters involving only 103 different modes. 
Moreover, 15 of these clusters have the simplest "triad" 
structure, formed by three modes. Importantly, there are 
only four isolated triads in the domain < m,i < 21, 
which is meteorological significant for the problem of cli- 
mate variability on intra-seasonal scale of about 10-100 
days (waves with £ > 21 have to short period to play 
significant role in this problem). 

The main physical message of our Letter is that so- 
called Intra-seasonal Oscillations (IOs) of the Earth at- 
mospheric flow can be rationalized as periodical energy 



exchange within the above mentioned four isolated triads 
of the planetary waves. IOs has been firstly have been 
detected in the study of time scries of tropical wind. 
Similar processes have been also discovered in the atmo- 
spheric angular momentum, atmospheric pressure, etc. 
Detailed analysis of the current state of the problem is 
presented in [10( and references therein; most part of the 
papers are devoted to the detection of these processes 
in some data sets, llj, [12[ and to the reproducing them 
in computer simulations with comprehensive numerical 
models of the atmosphere 
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Nevertheless, many as- 
pects of the IOs remain unexplained, e.g. the reason of 
IOs in the North Hemisphere is supposed to be topog- 
raphy, see e.g. [3] and no reason is given for IOs in the 
South Hemisphere, there is no known way to predict the 
appearance of IOs, etc. 

Our model considers IOs as intrinsic atmospheric phe- 
nomenon, related to a system of resonantly interacting 
triads of planetary waves, which is an example of wave 
mesoscopic system. The model is equally applied to the 
North and South Hemispheres, is independent (in the 
leading order) of the Earth topography, naturally have 
the period of desired order and allows to interpret the 
main observable features of IOs (see Sec. 5). 
1. Atmospheric planetary waves. These waves are 
classically studied in the frame of barotropic vorticity 
equation on a sphere: 



dAtp/dt + 2 dip/dX + J(V, AV>) = 



(1) 



Here ip is the dimcnsionless stream-function; velocity v = 
51 R [z x V-0] , with being the angular velocity of the 
Earth and z - the vertical unit vector; the variables t, <p 
and A denote dimcnsionless time (in the units of l/f2), 
latitude (— 7r/2 < <p < n/2) and longitude (0 < A < 
2tt) respectively; and J(a, b) are spherical Laplacian 
and Jacobian operators. The linear part of this equation 
has solutions in the form A 3 -.Y t 1 (A, (p) exp(iu)jt), u)j = 
—2mj/£j(£j + l). Integer parameters £j and (ij—mj) > 
are longitudinal and latitudinal wave-numbers of j-mode, 
they are equal to the number of zeros of the spherical 
function along the longitude and latitude. 

Assuming small level of nonlinearity, \Aj\ -C 1, we 
restrict ourselves by resonant interactions only. Under 
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FIG. 1: Four isolated triads are shown below the diagonal 
dash line; I - along horizontal axis and m < £ along the 
vertical one. The three clusters of two connected triads are 
shown above diagonal dash line (now £ - along the vertical axis 
and m < £ - along the horizontal one). The (meteorological 
significant) spectral domain is restricted by < m, £ < 21. 
Localized modes are shown by full circles, delocalized ones 
(with relatively small interaction amplitude) by empty circles. 



the resonance conditions for three modes: oj\ + L02 = 
uo 3 , mi + m 2 = m 3 , in which \£x — £ 2 \ < £3 < £\ + £2, 
and £\ + £2 + £3 is odd, the triad amplitudes Aj(i) varies 
in time according to equations 0]: 



NidAi/dt 
N 3 dA 3 /dt 



2ZN 32 A 3 A* 2 , N 2 dA 2 /dt = 2ZN 13 AlA 3 , 
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Table 1. For each triad the following data are given: all 
resonantly interacting modes, interaction coefficient Z, 
interaction latitude ifio (in grad), magnitude of the elliptic 
integral K([i), corresponding to the ESMRW December 
1989 data for 500 hPa initial energy distribution in a triad, 
and initial dimensionless energy Eo * 10 6 of each triad and 
resulting To values (in days). 



Here Nj = £j{£j + 1) , Nij = iV,- — Nj, and interaction co- 
efficient Z is an explicit function of wave numbers. This 
system conserves energy E and enstrophy H: 

E = E 1 +E 2 +E 3 , H = N\E\ + N2E2 + N 3 E 3 , (3) 

where the energy of j-modc is Ej = Nj\Aj\ 2 . 

2. Classification of the triads. Consider the struc- 
ture and properties of interacting resonant triads in the 
meteorological significant domain < m, £ < 21, where 
we found four isolated triads, denoted as Ai,... , A4, 
three "butterflies", i.e. clusters of two triads (denoted as 
txJi, M 2 and CXI3) that are connected by a common mode, 
and one cluster of 6 connected triads denotes as Kl. The 
structure of all isolated resonant triads and "butterflies" 
clusters is shown in Fig. [TJ Main information about the 
triads in the chosen spectral domain is given in left 4 
columns on Table 1: notations of the triads, three pair 



of J 



for each triad, the value of the interaction coef- 



ficient Z and the so-called "interaction latitude" ipo in- 
troduced in Ref. Columns 5-7 contain data which is 
necessary to compute period of resonant interactions and 
will be commented on further. 

We can interpret the latitude ipo as follows. The 
overlap of three wave-function in a triad, Z(X, ip) = 
r™ 1 (A,^)y™ 2 (A,^)F™ 3 (A,( < 3) shows a contribution to 
the interaction coefficient Z cx J 2(X, ip) dXdip from a 
particular location on the sphere. The overlap Z(X, ip) 



has a maximum at a particular latitude <po and a 
narrow latitudinal belt around ipo gives the main contri- 
bution to the global interaction amplitude Z. That is 
why ipo can be understood as the interaction latitude. 

3. General solution of the triad equations. Lin- 
ear change of variables Bi = onAi with a, being explicit 
functions on Nj allows to rewrite Sys. ((2|) as 

7£>i = 2ZB 2 B 3 , B2 = 2ZB*B 3 , B 3 = 2ZB\B2 ■ 
This system has two independent conservation laws 



h = 
h = 



\B 2 \ 2 



\B 3 \ 2 



(EN 1 ^H)N 23 /N 1 N 2 N 3 , (4a) 
|#i| 2 + \B 3 \ 2 = {EN 2 - H)N 13 /N!N2N 3 (4b) 



which are linear combinations of the energy E and en- 
strophy H . Direct calculations show that the general 
solution for Bi is expressed in Jacobian elliptic func- 
tions Bi = Bi, cn(T - r ), B 2 = B 2: oSn(r - r ), B 3 = 
B 3 odn (t — To), where Bj^ , tq , are defined by initial 
conditions and r = t/2 Z\Jl\l2 . Functions cn(r), sn(r) 
and dn(r) are periodic with the period AK (fi) , AK (fi) and 
2K(/j.) correspondingly, were K(/j,) 



d6 



/j 2 = min 



[i sin 



IV/iJ 



< 1 



Fig [5] illustrates the typical time dependence of all three 
dimensionless amplitudes of the triad Ai. One sees that 
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K{p) is a smooth function that changes slowly enough 
such for the wide region of the initial conditions it can 
be roughly considered as a constant. 
4. Period of triad oscillations. The period of energy 
exchange (measured in days) in the triads is given by 

T = 7r K(fi) I Z I 2 , that can be written a product of 

functions E and the ratio of the enstrophy to the energy 
h = H/E : 



T 



ME) 



T(E,h)=T (E)K((,)f(h), (5a) 

(5b) 



N 1 N 2 N 3 



2Zs/2E V N 21 ^/N 31 N 23 



f(h) = y/N 12 /2y/(Nr- h)(h - N 2 ) 



(5c) 



Here K depends on fx and, in its turn, fx depends on h 



H 2 (h) 



(h-N 2 )N 23 Ni — h) (N[ 
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Equa- 



(iVi-/i)7V 3 i {h-N 2 )N 23 . 
tion ([3]) show that possible values of h lie inside one of 
the two intervals N 2 > h > Ni or Ni > h > N 2 . With- 
out loss of generality we set N 2 > h > Ni, then maximal 
possible value h = N 2 is realized if Ei = £3 = 0, i.e. only 
the second mode is excited. The minimal value h = Ni 
is possible if E 2 = E 3 = 0, i.e. only the first mode is 
excited. In both cases, according to basic Eqs. ^ there 
is no time evolution, i.e. Ej = const, for j = 1, 2 and 3. 



This is in agreement with Eq. (|5c|) . according to which 
f(h) -► 00 for h -> N 2 or h -► N x . 

Function f(h), Eq. (|5c[) . has a minimum (equal to 
one) just in the middle of the interval N 2 > h > Ni 
at h = h + = (Ni + N 2 )/2. For this value of h 

{N ^ N3) iN *~ N ^ For isolated 



[i 2 (h + ) = min | 



(JVa-JVx) ' (JVa-iVa). 
triads of interest Ai , A2 , A3 and A4 , the values of 
p 2 {h + ) arc 0.93, 0.41, 0.97 and 0.45 respectively, with 
K(p) equal to 1.96, 1.27, 2.22 and 1.30. 

The less trivial case of infinite period corresponds to 
(j, = 1, which is realized at h = ft C rit = N\ + N 2 — N 3 . 
In this case B 3 (t) = tanhr and for r — * 00, S3 — ► 1 
and Bi and B 2 exponentially fast go to zero, i.e. for the 
specific value h = h cr n, the high mode B 3 exponentially 
fast takes energy from two low modes. This is possible 
only for three particular values of h: h — N\, N 2 and 
Ni + N 2 — N 3 . The h dependence of the period of the 
triads Ai , . . . , A4, is presented in Fig. [21 One sees, that 
the regions, where the period exceeds twice the minimal 
possible are very narrow, just few percent of the available 
interval of h. This means, that though theoretically for 
each triad we can always choose initial conditions in such 
a way that period will be large and even tend to infinity, 
the probability of this is very small. 

Indeed, for qualitative analysis we can think that in 
the turbulent atmosphere the probability to get some en- 
ergy from global disturbances to a particular planetary 
wave is independent from the state of other waves and 
this probability is more or less the same for each wave 



in a triad. If so, the probability V{h) to have initial 
conditions with some value of h has to be a smooth func- 
tion of h in the whole available interval N 2 > h > N\. 
Roughly speaking, we can approximate V{h) as the con- 
stant: V{h) ~ l/(N 2 - iVi), N 2 > h > iVi. With this 
approximation we can, for example, for triads A1-A4 es- 
timate the probability to have the period, twice exceeding 
the minimal one Tq(E), Eq. (|5bj) . as few percents. More- 
over, as one sees in Fig. 2, the typical value of the period 
T is about (1.4 ± 0.1)T for the triads Ai, A 3 and about 
(1.2 ± 0.1 )Tq for the triads A2, A4. This conclusion is 
in a qualitative agrement with the ECMWF (European 
Center for Medium-Range Weather Forcast) winter data, 
shown in Table 1, column 7. 

5. Intra-seasonal Oscillations as Resonant Triads. 

Our interpretation of IOs as dynamical behavior of 
Ai , A2 , A3 and A4 triads allows one to answer some 
questions appearing from meteorological observations 

• What is the cause of IOs in South Hemisphere ? The 
basic fact of our model is the existence of global nonlin- 
ear interactions among planetary waves, independent of 
the topography. 

• Why the period of so-called "topographic" oscillation 
in North Hemisphere is given as J^0 days by some re- 
searchers and 20-30 days - by other researchers? The 
variations in the magnitudes of the period are caused by 
different initial energy and/or initial energy distribution 
among the modes of the same triad. 

• How do the tropical and mid-latitude oscillations inter- 
act? Two mechanisms are possible: i) Triads with sub- 
stantially different interaction latitudes belonging to the 
same group, for instance, triads [(1,6) (2,14) (3,9)] and 
[(3,9) (8,20) (11,14)] of M exchange their energies through 
other modes of this group and belong correspondingly to 
the tropical and extra-tropical latitudinal belts, and ii) 
Isolated triads can interact via some special modes called 
active near-resonant modes Q. These modes have small- 
est resonance width with a given triad, and are them- 
selves parts of some other resonant triad. For instance, 
the mode (13,19) is a near resonant for A4 (with reso- 
nance discrepancy 8 = 0.16) and is resonant for A3. 

• Why do the intra-seasonal oscillations are better ob- 
servable in winter data? In summer, modes have higher 
energies, periods of the triads become smaller, and reso- 
nances with big enough resonance width can destroy the 
clusters. 

• How to predict these recurrent features? Amplitudes 
of the spherical harmonics with wave numbers taken 
fr om Table 1 have to be correla ted: (Ai(t)A 2 (t)A 3 (t)) — 
v /(|A 1 (t)|2)(|A 2 (i)| 2 )(|A 3 (t)| 2 ), sec Fig. 2. Magnitudes 
of the expected periods can be computed beforehand by 
the given explicit formulae. 

6. Conclusions 

• Our simple model provides main robust features 
of IOs in terms of resonance clusters consisting of three 
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t *. *. 

FIG. 2: Left panel: Time dependence of Ax, A2 and A3 [denoted by green solid, red dashed and blue dot-dashed lines 
accordingly] of the triad Ai, with n — 0.89, corresponding to observed data. Middle and Right panels: Dependence of the triad 
periods T(E,h)/To(E), Eq. (|5a[) . Middle panel: blue-dashed line corresponds to Ai triad, green-solid line - A3 triad. Right 
panel: Red dashed line - A2 and magenta solid line - A4 triad. Black dashed lines denote minimal period: T « 1.3Tb for Ai, 
A 3 and T « 1.18T for A 2 , A 4 . 



modes of atmospheric waves. 

• Energy behavior within the bigger clusters should be 
a subject of a special detailed study. Knowledge of clus- 
ter structure allows to simplify drastically their analysis. 
For instance, for "butterfly" cluster at least 6 real inte- 
grals of motions can be easily found. Universal method 
to construct isolated clusters and write out explicitly cor- 
responding dynamical equations for a wide class of meso- 
scopic systems is given in [l5j . 

• Our approach is quite general and can be used 
for studying many other mesoscopic systems, provided 
that explicit form of dispersion function 0(fc) is known 
(here k is the wave vector of plane systems with periodi- 
cal boundary conditions, or another set of eigenvalues in 
more complicated cases, like m, I for the sphere). Prop- 
erties of a specific mesoscopic system will depend on the 
1) form of 0(fe), 2) dimension of k, 3) number of conser- 
vation laws, 4) initial magnitudes of the conserved values 
(energy, enstrophy, etc.) and their initial distribution 
among the modes in the cluster. 
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